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Muon decay spectrum: leading logarithmic approximation 
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O (a 2 In 2 (m M /m e )) QED corrections to the electron spectrum and angular distribution in muon 
decay are evaluated. Impact on the determination of Michel parameters is estimated. Current 
theoretical uncertainty in the muon decay distributions is discussed. 
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I. INTRODUCTION 

The decay /i + — > e + v e v^ has been very accurately studied. Its total rate determines the Fermi constant G M 
describing the strength of the weak interactions jy. Standard Model predicts other features of this decay, such as 
the positron energy and angular distribution, positron polarization, and a variety of correlations between spins and 
momenta of the muon and its decay products , ||, ||[ |5| . Precise determinations of those observables test the Standard 
Model and can be used to search for possible "new physics" effects || . This motivates ongoing studies of the positron 
transverse Rl a nd longitudinal polarization || and angular /energy distributions J|, [To). Muon lifetime is also being 
remeasured |lTl [l2) . 

Current experiments are so precise that theoretical predictions must include radiative corrections beyond the first 
order in the fine structure constant a ~ 1/137.036. For the muon lifetime they are known []l3[ |l4|, |l5|, |l6|, [ij], [l^, |Tsfl 
through O (a 2 )- On the other hand, little is known about the positron energy distribution beyond the O (a) effects 
§§■ 

The experiment TWIST at Canada's National Laboratory TRIUMF is designed to measure the positron spectrum 
from polarized /i + decays with a precision of 10~ 4 |], [u|. To match this, and help search for "new physics" effects, 
the Standard Model prediction must include at least the leading O (a 2 ) effects in the positron distribution. 

We have examined the dominant effects at this order, which arise due to emission of collinear photons and e + e~ 
pairs. Such effects do not significantly affect the total muon lifetime, but they rearrange the positron spectrum. They 
are enhanced by two powers of the muon and electron mass ratio logarithm L = \n{rrvi/m\) ~ 10.66 (this roughly 
determines their order of magnitude, (a 2 /tt 2 )L 2 ~ 6 • 10~ 4 ). 

We are interested in the energy and angular distribution of positrons produced in the /i + decay. We normalize it so 
that it coincides with the differential width of the decay in the two lowest orders in a. A difference occurs at O (a 2 ) 
where an additional positron can appear due to pair production. The differential distribution of positrons (summed 
over e + spin states) in the polarized muon decay is 



- 1 j 3=r - r |F W - m »p,o(,)] . r„ s ^ (i + - -f- 
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where m e and m M are the electron and muon masses; 9 is the angle between the positron flight direction and the muon 
spin (for the decay, the sign of the cos 9 term should be reversed); P M is the degree of the muon polarization; E is 
the positron energy. 

Functions F(z) and G(z) describe the isotropic and anisotropic parts of the positron spectrum. They can be 
expanded in series in a, 

F(z) = f Boin (z) + |-A(z) + (|-) 2 h(z) + O (a 3 ) , (2) 
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and similarly for G(z). 

To match the precision of TWIST, the electron mass should be included at the Born level, 



/Born(z) = 6(1 + v'f VZ 
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where p, 77, £, and S are the so-called Michel parameters 



22I p3|] which depend on the Lorentz structure of the 



interaction responsible for the decay. In this paper we assume that the decay is caused by the Standard Model V — A 



interaction, for which p — 3/4, rj = 0, £ = 1, and <5 = 3/4. These values agree with present experimental fits 24 

p = 0.7518 ± 0.0026, 

T) = -0.007 ±0.013, 

= 1.0027 ± 0.0079 ±0.0030, 

5 = 0.7486 ± 0.0026 ±0.0028. 

In the massless limit (m e —* 0) we have 

/BornO) -> fo(z) = Z 2 (3 - 2z), gBorn («) ~> 50 [z) = Z 2 (1 - 2z) . 

Functions fi jl(J and 51 |^(| are also known with full dependence on the electron mass. However, at present it is 
sufficient to use their massless limit given in the Appendix. 

The O (a 2 ) effects are not yet known. They can be divided up into three parts according to the power of L, 
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and similarly for §2- In this paper we evaluate the leading logarithmic (LL) corrections /j/ L and g\ L . We divide them 
into contributions of pure photon emissions and of diagrams with e + e~ pairs, 



LL 
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In diagrams with e + e~ pairs we have to clarify the meaning of the variable z, whether it describes the energy of 
the "primary" positron or the one from the pair. (We can neglect interference in the LL approximation.) These two 
possibilities give rise to the so-called non-singlet (NS) and singlet (S) parts of pair corrections, 

,LL(e+e _ ) _ 2 LL(e + „LL(e+0 T.T. 2 

J 2 — g J 2 NS ± J 2 S ' 

Ingredients needed to evaluate the full LL effect in eq. (0) are given below in eqs. (]lj, |l5|) (photonic corrections) and 
m M) (pairs). 
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II. LEADING LOGARITHMIC APPROXIMATION 

The LL corrections can be found by convoluting the tree level spectrum (^[) with the positron structure function (SF), 
a solution of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution equations for QED. Analytical expressions for 
T>(w, (3) are known |25|, ^6], |?J up to the fifth order in a so that terms 0(a n L n ) can be found for n = 1, . . . , 5 (in this 
paper we treat n < 2). 

To find the various corrections outlined in the Introduction, we divide the SF into three parts: pure photonic, and 
non-singlet (NS) and singlet (S) e + e~ pair contributions, 

L — ' n! 

n— 1 

V s e+e .(w,l3) = ^R(w)+0(a 3 L% 0=£-(L-l). (9) 
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The components Z? 7 and Z> + _ correspond to those Feynman diagrams in which the registered positron belongs to 
the same fermionic line as the initial one. If the registered positron arises from a pair production, its LL contribution 
to the energy spectrum is described by the singlet function T>^ + , Functions relevant for our work are 
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Higher order expressions and further details on the SF formalism can be found in p||, ^7], g8|, ^9|, |3(J . 
To find the LL corrections we use convolution, defined by 

11 l 

A(»)®B(z) = J dw J dw' S(z - ww')A{w)B{w') = J ^A(w)B\J^J. 

z 

For example, to reproduce the first order LL correction, we convolute the Born-level spectrum with P^ 

A LL (z) = pV(*)®f (z) = ~+2z-4z 2 +?-z 3 + 2z 2 (3-2 Z )ln—, 

6 3 z 

g\\z) = - \-Az 2 + \z* + 2z 2 (l-2z)lni^. 

6 3 z 



(11) 

(12) 
(13) 



These formulas coincide with the LL parts of the full O (a) results given in the Appendix. A comparison of the LL 
and full first order functions is presented in Fig. |l|. We see that the LL approximation gives the bulk of the O (a) 
correction, especially in the region of intermediate and large values of z, relevant for TWIST. 
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FIG. 1: Values of functions /i and gi versus z: exact results (eq. (Al)) and LL approximations (eqs. ( |12[|l.'j| )). 
For the second order photonic LL corrections, we convolute with P^> , 
f^\z) = P^{.)®fo{z) 



4z 2 (3 - 2z)$(z) + 



— + 8z-16z 2 + — z 3 )ln(l-z) 
3 3 



32 



6 



11 17 8 , 32 



—x - 2z + 8z 2 - — z A lnz + — + — z + -z z - — z , 



36 6 



(14) 



5 2 L(7; (z) = 4z 2 (l - 2z)$(z) + I -- - 16z 2 + — z A ) ln(l - z) 



.2 , 32 3 



4 



32 



5z — z 



7 7 8 2 32 

36 -6 z+ r-Y z 



$(z) = Li 2 



1 - z 



In" 



1 — 2 7T 



G 



Li 2 (a:) 



dy 



In(l-y). 



(15) 



In the same manner we can get the third order photonic contributions. 
Integrals of the LL photonic contributions vanish, 

11 11 

dz A LL (z) = / dz f^\z) = fdz gp(z) = fdz g^\z) = 0, 



(16) 



in accord with the theorem about the cancellation of mass singularities ]3l], [3^] . 

A numerical illustration of our results for the relative size of the second order LL photonic corrections is given in 
Fig. where we plot the relative correction defined as 
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FIG. 2: Relative size of the second order LL photonic corrections for as a function of z, plotted for three values of cos 8. 

In O (a 2 ) the positron distribution has a contribution from real and virtual e + e~ pairs. Virtual effects of heavier 
fermions are negligible |37j |. Pair correcttions are found by convoluting fo and go with (non-singlet) and R 
(singlet), 



/ 2 L N5 + °« = A LL W> g\\f e ~\z)=g^{z) (see eqs. (12,y 
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(19) 



In Table I, numerical results are presented for the quantity e \ which gives the relative size of the second 

order LL pair correction with respect to the Born distribution (for P M = 1), 



^llo+o _ f dw fp(z/w) - cos 9g (z/w) 
w f (z) - cos 9g {z) 



Wmin = max < z, 



z + y 



(20) 



where y is the cut on the maximal energy fraction of the real pair, -E pa i r < ym^/2. Both the singlet and the non- 
singlet pair contributions are taken into account. It means that we simulated the situation, where an observation of 
two positrons is treated as a pair of simultaneous /i + decays. 
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TABLE I: Leading logarithmic pair correction lO 4 ^ 1 ^ e ' for cos (9 = — 1. 



z\y OA 02 03 04 0.5 

0.05 5.60 20.62 51.01 99.44 164.88 

0.1 1.90 5.15 10.03 16.83 25.32 

0.2 0.46 1.66 2.83 4.07 5.35 

0.3 -0.10 0.73 1.36 1.91 2.36 

0.5 -0.72 -0.14 0.17 0.35 0.40 

0.7 -1.18 -0.80 -0.71 -0.71 -0.71 

0.9 -2.01 -2.01 -2.01 -2.01 -2.01 



In general, effects due to e + e~ pairs depend on experimental conditions and cuts for events with several charged 
particles in the final state. 

If one is interested in the total LL effect due to real and virtual e + e~ pairs to the muon decay width, one has to 

drop the singlet pair contributions f^s^ & ( z ) an d #2 & ^(z), to avoid double counting of real pairs. Cancellation 
of the leading logarithms in the total decay width is guaranteed by ([l6]) . 



III. EXPONENTIATION 



Toward the energy spectrum end point {z — ► 1), the first order correction and diverge. This phenomenon, 
discussed in Ref s. lira , |33| , is a signal to look beyond the first order approximation. One can use the Yennie-Frautschi- 
Suura theorem ]34j to re-sum the divergent terms and convert them into an exponential function. Exponentiated 



representations of the SF can be employed to re-sum parts of the leading logs to all orders in aL 28, |35[. The 



exponentiation for the muon decay has been criticized ]36| , because the large logarithmic terms contain a mass 
singularity: not all large logs disappear in the terms supplied by the exponent after an integration over the energy 
like in eq. (^6|). 

As has been discussed in Q , the validity of exponentiation is limited to the region near the end of the spectrum. 
One can see from the O (a) results that the exponentiation can be relevant only in a very small range where z differs 
from 1 by about 10 -10 (the correction is about —50% at 1 — z = 1CP 10 ), which is much less than the experimental 
resolution. For this reason we leave our results in the un-exponentiated, fixed order form. 

The end region of the spectrum is usually excluded from the fits of Michel parameters. This is done to reduce 
the uncertainty due to the finite energy resolution, which is most important in this region. Indeed, the shape of an 
experimentally observed spectrum is a convolution of the "true" spectrum with a resolution function. In the "bulk" 
part the effect of the finite resolution is very small. However near a sharp edge (with the width much less than the 
width of the resolution function) the shape of the convoluted spectrum is defined mainly by the resolution function 
and not by the spectrum. Exclusion of the end point region helps also to reduce the theoretical uncertainty because 
this is where the unknown higher-order corrections are expected to be the largest. 



IV. CONCLUSIONS 



To estimate the effect of the second order correction on values of Michel parameters measured in an experiment, 
we generated a 2D distribution in z and cos 9 according to the RC-corrected spectrum and fitted with a spectrum 
without the corrections. 

10 9 toy Monte Carlo "events" {z,cos0} were produced by sampling the 2-dimensional spectrum ([[]) taking into 
account the complete first order corrections and the second order LL photonic ones. This level of statistics is expected 
to be accumulated in the TWIST experiment [^0|. The acceptance-rejection method and the Mersenne Twistor 
random number generator Q was used. Events passing the "acceptance cuts" 0.34 < |cos#| < 0.98, 0.4 < z and 
the cut z < z max were filled into a 2D histogram. z max varied between 0.96-0.995. The cuts roughly represented 
acceptance of the TWIST detector |l(J . The region of z close to 1 was excluded to avoid the issue of the experimental 
resolution, discussed above. Finally, a maximum loglikelihood fit of the spectrum without the second order RC to the 
histogram was done, p, rj, £, S, and the global normalization were the 5 free parameters of the fit. We put = 1 
both in the generation and in the fits. Binning of the histogram was chosen sufficiently fine, so that repeating the 
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procedure with several times smaller bins gave the same results. Self-consistency of the method was checked by fitting 
the histogram with the full spectrum. Original values of Michel parameters were reproduced within the errors. 

We have observed statistically significant deviations of the fitted Michel parameters from their original values when 
doing the fits without the and g\ L ^ contributions. That emphasizes the importance of using a precise enough 

theoretical spectrum shape for extracting values of Michel parameters in an experiment. 

Shifts of Michel parameters due to radiative corrections depend on the fit region, and perhaps on other factors not 
considered here. For example, one may want to take into account the effect of the finite experimental resolution in the 
bulk part of the spectrum. Fig. |^ demonstrates dependence of the shifts on the upper energy limit of the fit region. 

For a realistic value of the cut, 2 max = 0.97, the shifts of Michel parameters due to the second order LL corrections 
are of the order 

Ap ~ 11 -lCT 4 , 
Ar) ~ 350 -lCT 4 , 
A£ ~ 3-10~ 4 , 

AS ~ 4-10" 4 . (21) 
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FIG. 3: Shifts of Michel parameters due to the second order radiative correction for different upper energy cuts. Horizontal 
axis is z m ax, vertical axis is the difference between a reconstructed parameter and its Standard Model value. The points are 
correlated since they are obtained from the same data set. 

The relatively large shift of rj is due to the fact that it enters the spectrum with a small coefficient zq ~ 10~ 2 . If 
we fix r} during the fits to its SM value 0, the shifts due to O (a 2 ) effects become 5 ■ 10~ 4 for p, —3 • 10~ 4 for £, and 
3 • 10 -4 for S. Comparing the results for free and fixed r\ we observe a strong correlation of r\ with p and £. 

Fig. |^ indicates a rather strong dependence of all Michel parameters on z max . It follows from the peaked behavior 
of the radiative corrections when z is close to 1. The fitting procedure tries to compensate the effects of radiative 
corrections and cuts by adjusting Michel parameters and the global normalization. 
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The remaining theoretical uncertainty in the muon decay spectrum, eq. (|l]), is due to unknown contributions 
O (a 2 L) and 0(a n L n ), n > 3. Non-logarithmic terms of the order (a/71") 2 ~ 5.4 • 10~ 6 are expected to be small 
(compared with the 10 -4 precision tag). The sub-leading contributions O (cn 2 L) are the main source of the remaining 
uncertainty. The magnitude of the corresponding photonic contribution can be estimated using the known second 
order photonic LL correction, 

5 ^~^ h ^^0.3 6^\ (22) 

where we multiplied by 3 to account for the unknown coefficient of sub-leading terms. The leading contributions from 
higher orders (n > 3) can be estimated using again the second order LL result, 



S^™ ~ 3j-LS^ (7> ~ 0.03 S^> . (23) 



<j; LL(7) _ T X Lh(-y) ^nnq 

e + e~ pair contributions are typically smaller than those of the photons, at least in the LL. 

These estimates of unknown radiative corrections can be converted into theoretical uncertainties a th of Michel 
parameters, of the order of a third of the shifts in eq. ( |2l| ) we found by including the LL photonic corrections. With 
z-m&x — 0.97 we find 

a th {p) = 3-f0~ 4 , 
a th (ri) = fOO -f0- 4 , 



th/jn _ 1 in-4 



<7 tn (<5) = l-MT*. (24) 

The conditions and the fitting procedure in a concrete experiment can be different from the ones described above. The 
actual size of the effect of radiative corrections on Michel parameters can be derived there in a similar way, starting 
from the analytical formulas for theoretical predictions and applying specific experimental conditions. 
The planned accuracy of the TWIST experiment j|, [l0| is 

a cxp (p) = 1-10" 4 , 
cr cxp (rj) = 30 -fO" 4 , 

<t cxp (6Pm) = 1.3 -to- 4 , 

a cxp ((5) = 1.4 -fO^ 4 . (25) 



Clearly, effects of the second order LL radiative corrections, eq. (|2l|), have to be taken into account at this level of 
accuracy. In order to further reduce the theoretical uncertainty, the next -to-leading second order corrections should 
be evaluated as well. Work on this is in progress. 
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APPENDIX A: FIRST ORDER CORRECTIONS 

O (a) corrections to the muon decay spectrum, without terms suppressed by m 2 /m 2 l , read [fl6[ 

f x (z) = (L-l)fl jL (z) + 2z 2 (3-2z)R 1 (z) + — - [(10 + 34z - 32z 2 ) lnz + 5 - 27z + 34z 2 l , (Al) 

6 

9l (z) = (L — l)g^(z) + 2z 2 (l - 2z)R 1 (z) - 1 + 2?Z " ~ ^ lnz - i-^ (7 - 13z - 30z 2 ) - lii-Zfl! m (i _ z ). 



3z 



Rxiz) = -2Li 2 (l-z) + lnzln(I-z)-21n 2 z-^ — - ~~ A 

z 4 
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fi h (z) and g\ h (z) are denned in eqs. (||JT|). 
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